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Abstract 

We consider the eigenvalue problem for the Laplace operator in a planar domain which can be decomposed 
into a bounded domain of arbitrary shape and elongated "branches" of variable cross-sectional profiles. When 
the eigenvalue is smaller than a prescribed threshold, the corresponding eigenf unction decays exponentially 
along each branch. We prove this behavior for Robin boundary condition and illustrate some related results 
by numerically computed eigenfunctions. 

1 Introduction 

The geometrical structure of Laplacian eigenfunctions has been thoroughly investigated (see the review [1] 
and references therein). When a domain can be seen as a union of two (or many) subdomains with narrow 
connections, some low-frequency eigenfunctions can be found localized (or trapped) in one subdomain and 
of small amplitude in other subdomains. Qualitatively, an eigenfunction cannot "squeeze" through a narrow 
connection when its typical wavelength is larger than the connection width. This qualitative picture has found 
many rigorous formulations for dumbbell shapes and classical and quantum waveguides ^2-nl2j. Numerical and 
experimental evidence for localization in irregularly-shaped domains was also reported p!3H2T] . 
In a recent paper, we considered the Laplacian eigenvalue problem 

Au -\- Xu = in 

(1) 

u = on dQ, 

for a large class of domains Q in {d = 2,3, ...) which can be decomposed in a "basic" bounded domain V 
and a branch Q of a variable cross-sectional profile [10 . We proved that if the eigenvalue A is smaller than the 
smallest eigenvalue fi among all cross-sections of the branch, then the associated eigenfunction u exponentially 
decays along that branch: 

Theorem 1. Let Q C (d = 2,3, ...y) be a bounded domain with a piecewise smooth boundary and let 
Q{z) = ft n {x e : xi = z} the cross-section of ft at xi = z e R by a hyperplane perpendicular to the 
coordinate axis x-\ (Fig. [7p. Let 

zi = mf{z e R : Q{z) ^ 0}, Z2 = sup{z G R : Q{z) / 0}, 

and we fix some zq such that Zi < zq < Z2. Let ii{z) be the first eigenvalue of the Laplace operator in Q{z), with 

Dirichlet boundary condition on dQ{z), and ji = inf l^{z). Let u be a Dirichlet-Laplacian eigenfunction in 

ze{zo,z2) 

Q satisfying and A the associate eigenvalue. If \ < ji, then 

\Ml2{Q{z)) < II^IIl2(Q(^o)) exp(-/3V/i- a {z - zq)) {z > zo), (2) 

with P = 1/a/2. Moreover, if (ei • n{x)) > for all x G dft with xi > zq, where ei is the unit vector (1, 0, 0) 
in the direction xi, and n{x) is the normal vector at x e dft directed outwards the domain, then the above 
inequality holds with j3 = 1. 
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Figure 1: Two examples of a bounded domain Q = with a branch of variable cross-sectional profile 

Q{z). When the eigenvalue A is smaller than the threshold /i, the associated eigenfunction exponentially decays 
in the branch Q2 and is thus mainly localized in Qi. Note that the branch itself may even be increasing. 



In this theorem, a domain Q is arbitrarily split into two subdomains, a "basic" domain Qi (with xi < zq) 
and a "branch" Q2 (with xi > 2:0), by the hyperplane at Xi = zq (the coordinate axis Xi can be replaced by any 
straight line). Under the condition A < /i, the eigenfunction u exponentially decays along the branch 0.2- Note 
that the choice of the splitting hyperplane (i.e., 2^0) determines the threshold fi. Since fi is independent of the 
basic domain V, one can impose any boundary condition on dQi (that still ensures the self-adjointness of the 
Laplace operator). In turn, the Dirichlet boundary condition on the boundary of the branch Q2 was relevant. 
Many numerical illustrations for this theorem were given in [10 . 

Remark 1. It is worth stressing that the sufficient condition A < /i involves purely spectral information: the 
eigenvalue A in the whole domain and the smallest eigenvalue fi over all cross- sections. A simple geometrical 
condition on the basic domain Qi can be formulated through the inradius p of Qi (or Q), i.e., the radius of the 
largest inscribed ball Bp. Since Bp C Vt, the first Dirichlet eigenvalue A is bounded as X < Xi{Bp) = , 

where j d_^ is the first positive zero of the Bess el function Jd_^{z). A sufficient geometrical condition for getting 
exponentially decaying eigenfunction is then 



For planar domains, the inequality yields p/h > Jo/tt where b is the largest width of the branch, and jo 2.4048. 
This inequality includes only the inradius of Qi (or Q) and the largest width of the branch, while the length of 
the branch can be varied arbitrarily. For instance, the localization in the basic domain Qi may hold even when 
the area of the branch 1^2 '^s arbitrarily large, as compared to the area of Vti. 

Remark 2. For higher dimensions (d > 3), the localization may sound even more striking, as the %ranch^' 
has to be ^^narrow" only in one direction (Fig. [7p. For instance, if the branch 0.2 has a constant width b in one 
direction, then the smallest eigenvalue p in its cross- sections is greater than 7r^/6^. If the inradius p of is 
greater than bjd_^/Ti then the inequality 151) holds, and at least the first eigenfunction is localized in l^i. In the 
three-dimensional space, ji = tt so that the inradius has to be just greater than b: p > b. A simple example 
is a domain decomposed into the unit cube Qi = {{x^y^z) gM^ : 0<x<l, 0<?/<l, 0<z<l} and 
a parallelepiped Q2 = {{x^y, z) e : 1 < x < L^, — Ly < y < Ly^ < z < b}. When b < 1/2, the first 
eigenfunction is localized in the cube fti, whatever the lateral spatial sizes and Ly of the %ranch" are. 

In the remainder of the paper, we extend the above result to the Laplace operator in planar domains with 
Robin boundary condition. We also provide several numerical illustrations of localized eigenfunctions in planar 
domains in Sec. [31 

2 Extension for Robin boundary condition 

We consider the eigenvalue problem for the Laplace operator in a planar domain ft with Robin boundary 
condition on a piecewise smooth boundary dQ: 



where is a nonnegative function, and d/dn is the normal derivative directed outwards the domain. In that 
follows, we prove the following 

Theorem 2. Let Q = Qi U Q2, where Qi cM? is a bounded domain, and 



(3) 



Au -\- Xu = in 1], 
Ou 

- — \- hu = on 917, 
on 



(4) 



^2 = {ix,y) G : < a; < a, yi{x) <y< y2ix)} , 



(5) 
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is a branch of length a > and of variable cross-sectional profile which is defined by two functions ^1,^2 ^ 
(7"'^([0,a]) such that y2{ci) = yi{ci), yi{x) > and y2{x) < 0. Let u and A be an eigenfunction and eigenvalue 
of ft satisfying the eigenvalue problem Uv, with a nonnegative function h. We define ji = inf where 

I I xo<x<a 

is the first eigenvalue of the Laplace operator in the cross- sectional interval [yi{x)^y2{x)]: 

"^"{y) + l^i{x)v{y) = {yi{x) <y < y2{x)) 

v\y) - hi{x)v{y) = {y = yi{x)) (6) 
v\y) + h2{x)v{y) = {y = y2{x)) 

where 

hi{x) = h{yi{x))^l + [y'^{x)]^ (i = l,2). (7) 

If \ < 11, then 

\ML^{^a{x)) < \\u\\L2mxo)) exp(-/3V/i- A {x - xq)) {x > xq), (8) 
where /3 = 1/V2 and ^{xq) = {{x^y) eM? : xq < x < a, yi{x) < y < y2{x)}. 

Proof The proof relies on Maslov's differential inequality and follows the scheme that we used in [10^ for 
Dirichlet boundary condition. We consider the squared L2-norm of the eigenfunction u in the "subbranch" 
ft{xo): 

a y2ix) 

Iixo)= I dxdy = I dx I dy uHx,y) 

Qixo) xo yi{x) 

and derive the inequality for its second derivative: 

r(xo)>2(/i-A)/o(xo). (9) 

(i) From the first derivative 



I\xo) = - J dyu^{xQ,y), 



y2(xo) 
yiixo) 

we obtain 

2/2(3^0) y2ixo) 



dy u— - y2{xoW{xo,y{xo)) ^ y[{xoW{xo,yi{xo)) > -2 / dy u — , 



yi{xo) yi{xo) 
where we used the conditions y2{x) < and y'lix) > 0. Taking into account that 

y2ixo) y2{xo) 

f 1 9^ f 1 9^ f ir. 9^ f m 1 2 

— ay u-^ = — ay u- — V I ab u- — h j ab an 

yi{xo) yi{xo) S{xo) S{xo) 

dxdy dYv{uVu) -\- dS hu^ 



^{xq) S{xo) 

j dxdyiVu,Vu)+ I dxdyuAu+ J dS hu^ 

ft{xo) ^(xo) S{xo) 

dxdy iVu,Vu) - X I dxdyu^^ j dS hu^ 



^{xq) ^{xq) S{xo) 



where 



5 = 5iU52, S, = {{x,y)e^^ : {)<x<a, y = y,{x)} (i = l,2) 
is the "lateral" boundary of 1^2, we obtain 
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"{xq) > 2 j dS hu^ ^2 j dxdy{^^ - 2A j dxdy u 

S{xq) n{xo) n{xo) 

a y2{x) 2 

= 2 j dx^2{x) u^{x,y2{x)) ^hi{x) u^{x,yi{x)) ^ j dy - A^i^ | 



xo yi{x) 
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where hi{x) is defined by Eq. 

According to the Rayleigh principle, the first eigenvalue of the eigenvalue problem ([6| on the cross- 

sectional interval [yi{x),y2{x)] can be written as 



= inf 



h2{x) u'^{x,y2{x))^hi{x) u'^{x,yi{x)) 



y2{x) 



yiix) 



ueH^i[yiix),y2ix)]) y2{x) 

J dyu^{x,y) 
yiix) 



(10) 



from which we get 



y2ix) a y2{x) 



I (^o) — ^ J dx{iii{x) — A) j dy u (x, y) > 2(/i — X) J dx J dy u (x, y) = 2(/i — A)/(xo), 

XQ y^{x) Xo y^{x) 

where /i = inf fii{x). That completes the first step. 

xo<x<a 

(ii) We easily check the following relations: 

/(a) = 0, /'(a) = 0, I{xo)^0 (0 < < a), /'(xq) < (0 < < a). (11) 

Note that the second relation relies on the assumption that yi{a) = ^2(<^)- 

(iii) From the inequality (|9| and relations (11), an elementary derivation implies the inequality (|8|. In fact. 



one multiplies ([9]) by /'(xq), integrates from to a, takes the square root and divides by I{xo) and integrates 
again from xq to x (see [10 for details). □ 

The statement of Theorem [2] for Robin boundary condition is weaker than that of Theorem [l] in several 
aspects: 

• Theorem [2] employes an explicit parameterization of the branch through smooth height functions yi and 
^2; in particular, the statement is limited to planar domains. 

• The branch has to be non-increasing (conditions yi{x) > and y2{x) < 0) and vanishing at the end 
(condition yi{a) = y2{o))' 

• The inequality ([8| characterizes the L2-norm of the eigenfunction in the distant part of the branch, Q.{x)^ 
while the inequality ([2| provided an estimate at the cross- section Q{x). 

• The decay rate in Eq. ^ involves the coefficient P = 1/V2 while the inequality ^ for non-increasing 
branches was proved for /3 = 1. 

These remarks suggest that the statement of theorem [2] can be further extended while certain conditions may 
be relaxed. 

We also note that the solution of the eigenvalue problem ([6| has an explicit form 

v{y) = Cism{ay) -\- C2 cos (ay), (12) 

with two constants ci,C2 and fii{x) = a^, while the boundary conditions at the endpoints y = yi{x) and 
y y2{x), 

ci[—hi sin(a?/i) + acos(a?/i)] + C2[—hi cos(a?/i) — asm.{ayi)] = 0, 
ci [h2 sm{ay2) + a cos{ay2)] + C2 [/12 cos{ay2) - a sin{ay2)] = 0, 

yield a closed equation on a: 

(a^ + hih2) sina{y2 - yi) + a{hi + /12) cosa{y2 - yi) = (13) 

(here hi ^2 and yi^2 depend on x). This equation has infinitely many solutions that can be found numerically. 
The first positive solution will determine /ii(x). 

3 Illustrations 

In order to illustrate the geometrical structure of Laplacian eigenfunctions, we compute them for several simple 
domains. For all considered examples, we impose Dirichlet boundary condition for the sake of simplicity. 
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Figure 2: First three Dirichlet Laplacian eigenfunctions for the planar domain Q with sine-shaped branches, 
with 6 = 1, L = 1.54 and a = 1, a = 5 and a = 20. The first eigenfunction exponentially decays along the 
branch while the second and third eigenfunctions do not. The localization occurs in spite of the fact that 
the area of Qi presents only 7.15% of the total area for the last domain with a = 20. Note that the third 
eigenfunction for a = 5 is localized at the end of the branch 1^2 • 

3.1 Sine-shaped branches 

We consider the planar domain Q composed of a basic domain Qi (square of side L) and a branch Q2 of constant 
profile: 

n = {(x, y)eR^ : xe (0, a), y e {f{x)J{x) + b)}. (14) 

For this example, we choose f{x) = sin(x), fix 6 = 1, L = 1.54 and take several values for the length a. Since 
the inradius of the square fti is greater than Jo/tt, the first eigenvalue A in these domains is smaller than /i = tt^ 
for any length a so that the first eigenfunction should be localized in and exponentially decay along the 
branch 1^2 • This behavior is illustrated on Fig. [2] 

3.2 Star-shaped domains 

Figure [3] shows the first five Dirichlet Laplacian eigenfunctions for a "star-shaped" domain which is formed by 
a disk with many elongated triangles. The inradius of this domain is much greater than the largest width of 
triangular branches that implies localization of the first eigenfunction. One can see that all the five eigenfunctions 
are localized in the disk and exponentially decay along the branches. 

3.3 Elongated polygons 

As we discussed at the beginning, the separation into a basic domain and a branch is conventional. We illustrate 
this point by showing the exponential decay of the first Dirichlet Laplacian eigenfunction in elongated polygons 
for which the ratio between the diameter and the inradius is large enough. We start by considering a right 
triangle then extend the construction to general elongated polygons. 

We consider a rectangle of sides a and b {a > b) on which a right triangle ft with legs c and d is constructed 
as shown on Fig. [4j Note that the triangle is uniquely defined by one leg (e.g., d), while the other leg is 
c = ad/{d — b). The vertical line ai x = a splits the triangle Q into two subdomains: Qi (a trapeze) and (a 
triangle). For fixed a and 6, we are searching for a sufficient condition on d under which the eigenfunction u 
satisfying the eigenvalue problem ([T]), is localized in fti and exponentially decays along the subdomain (^2- 
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Figure 3: Localization of the first five Dirichlet Laplacian eigenfunctions in a domain Q with 51 branches. 




Figure 4: Left: A right triangle Q is decomposed into a trapeze Qi and a right triangle 1^2- Right: An 
elongated polygon P with n vertices {0, 5, A, Pi, P2, . . . , P/c}, where k = n — 3. Here, Qs is the polygon 
including n — 1 vertices {0, A, Pi, P2, . . . , Pfc}. 
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Lemma 1. Let Q be a right triangle Q defined by fixed a, 6, d > (a > b) (Fig. Let ( = (d/b) — 1, Pa{0 ci^d 
Pb{C) be two explicit polynomials defined by Eq. (16), and Co ^ 0.0131 is the zero of Pa{C)- If both inequalities 



a 



2 



OCo, (15) 

are fulfilled, then the first eigenfunction u of the Laplace operator in with Dirichlet boundary condition 
exponentially decays in 0.2- 

Proof. The proof of the exponential decay reUes on Theorem [T] For the "branch" (^2, the largest width is b so 
that the threshold ji = 7r^/6^. Our goal is therefore to find a sufficient geometrical condition to ensure that the 
first eigenvalue Ai is smaller than ji. This condition can be replaced by a weaker condition 71 < 7r^/6^ for the 
first eigenvalue 71 of the Laplace operator in the trapeze Vti with Dirichlet boundary condition. 
The eigenvalue 71 can be found from the Rayleigh's principle as 

71 = mf 7(v), 7(v) = — . 

Taking a trial function 

v{x,y) = yfy - d^ - — -x] sin(7rx/a), 



which satisfies Dirichlet boundary condition on the boundary of l^i, we look for such conditions that 7(1') < 
TT^b^, i.e., 



= -p-(^,^)L2(Qi) - (Vv, Vv)l2(Qi) > 0. 
The direct integration yields 

where hi = a/6, and Pa(C), Pb{C) are two polynomials of the fifth degree: 

5 5 



with the explicit coefficients: 



As 


= 27r4- 


157r2 + 


45 « 91.7741 




= 2-K^ + 


157r2 - 45 « 297.8622 


A4 


= 6(27r4 


-IOtt^ 


+ 15) « 666.7328 


B4 


= 6(27r4 


+ IOtt^ - 15) « 1671.0854 


A3 


= 30 (tt^ 


-47r2^ 


- 3) « 1827.9202 


B3 


= 30 (tt^ 


+ 37r2) Si 3810.5371 


A2 


= 20(27r 




+ 9) « 2299.8348 


B2 


= 20(27r 


* + 37r2) Si 4488.5399 


Ai 


= 30 (tt^ 


-67r2) 


w 1145.7439 


Bi 


= 307r4 f 


« 2922.2727 


Ao 


= 12(7r4 


-IOtt^ 


) « -15.4434 


Bo 


= 127r4 f 


« 1168.9091 



Note that all Bj > and Aj > except for Aq < 0. From the fact that Aj < Bj, one has (5(C) < for all C > 
when = 1 (i.e., a = b). We have therefore two parameters, C and a^, which determine the sign of Q and thus 
the exponential decay. Since Pb{C) > for all C > 0^ the condition (5(C) > is equivalent to two inequalities: 



.2 . Pb{0 



P^(C)>0, « >;^- (17) 

One can check that Pa (Co) = at Co ^ 0.0131 and Pa{C) > if and only if C > Co that completes the proof. □ 

We remind that this condition is not necessary (as we deal with an estimate for the first eigenvalue). For 
given a and b (i.e., a^), the above inequalities determine the values of C (and thus the leg d) for which localization 
occurs. Alternatively, one can express a and b through the legs c and d (and parameter C) as 

cC , d 

a = T — Ti = 



C + 1' C + 1' 

from which k = c(/d. For given c and d, one can vary ( to get a family of inclosed rectangles (of sides a and 
b). The above inequahties can be reformulated as 

C>Co ^ b< 



Co + 1' 
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Q>0 



Q<0 



2 4 6 8 10 

e 

Figure 5: The diagram in the space of parameters ( and for positive and negative signs of Q which correspond 
to locaUzation and non-locahzation regions. For a given (e.g., n = 2 shown by horizontal dotted hne), one 
can determine the values of for which localization occurs. 
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Figure 6: Several Dirichlet Laplacian eigenfunctions in (a) the right triangle with a = 2, 6 = 1 and ( ^ 0.32 
for which c = 8.25 and d = 1.32; (b) the right triangle with a = 4, 6 = 1 and ( ~ 0.08 for which c = 61.14 
and d = 1.07; and (c) elongated hexagon. In all these cases, the first eigenvalue Ai is smaller than tt^, while 
the associated eigenfunction decays exponentially along the "branch" (^2- The other eigenfunctions are also 
concentrated in l^i. 



The function f{Q can be checked to be monotonously decreasing so that the last inequality yields 

Of-Hc/d) ^ ^<pT(4^, (18) 

where denotes the inverse of the function /(C)- This condition determines the choice of the inscribed 
rectangle (the size b) for a given triangle. 

For the "worst" case c = d^ for which a numerical computation yields /~^(1) ^ 1.515, one gets 

This example shows that one can always inscribe a rectangle in such a way that A < ix^ j}? . However, the 
"branch" Q.^ in which an exponential decay of the eigenfunction is expected, may be small. Figure [6] illustrates 
these results. 

Remark 3. Any enlargement of the subdomain Qi on Fig. ^further diminishes the eigenvalue 71 and thus 
favors the exponential decay in ^2- In particular, for each positive integer n (n > 3), one can construct elongated 
polygons of n vertices for which the first Dirichlet Laplacian eigenfunction is localized in (Fig. Figure 
\Ejp shows first eigenfunctions in elongated hexagons. 
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